Quantum tunneling induced Kondo effect in single molecular magnets 



o 
o 



C. Romeike, M. R. Wegewijs, W. Hofstetter, and H. Schoeller 
Institut fur Theoretische Physik A, RWTH Aachen, 52056 Aachen, Germany 

We consider transport through a single-molecule magnet strongly coupled to metallic electrodes. 
We demonstrate that for half-integer spin of the molecule electron- and spin-tunneling cooperate 
to produce both quantum tunneling of the magnetic moment and a Kondo effect in the linear 
conductance. The Kondo temperature depends sensitively on the ratio of the transverse and easy- 
axis anisotropies in a non-monotonic way. The magnetic symmetry of the transverse anisotropy 
imposes a selection rule on the total spin for the occurrence of the Kondo effect which deviates from 
the usual even-odd alternation. 

PACS numbers: 72.10.Fk, 75.10.Jm, 75.30.Gw, 75.60.Jk 
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Introduction. Single-molecule magnets (SMMs) such 
as Mni2 or Fes have been the focus of intense experimen- 
tal and theoretical investigation [lj . These molecules are 
characterized by a large spin (S > 1/2), easy-axis and 
transverse anisotropies, and weak intermolecular inter- 
action. Molecular-crystal properties are due to an en- 
semble of single molecules and exhibit quantum tunnel- 
ing of magnetization (QTM) on a mesoscopic scale. Re- 
cently, a single molecule magnet (Mn.12) was trapped in a 
nanogap [2j,|3( and fingerprints of the molecular spin were 
observed in electron transport. Furthermore, transport 
fingerprints of QTM were predicted Q| when the indi- 
vidual excitations can be resolved by the temperature. 
Using easy-axis anisotropy for magnetic device operation 
was also proposed H . These works focused on the regime 
where single electrons charge and discharge the molecule 
through weak tunneling. 

In this Letter we investigate linear transport through a 
half-integer spin SMM deep inside the blockade regime |fj 
where the charge on the molecule remains fixed. A strong 
tunnel-coupling to the metallic electrodes induces spin 
fluctuations and allows the magnetic moment to tun- 
nel. This is remarkable, since for an isolated SMM 
with half-integer S this is forbidden by time-reversal 
(TR) symmetry. At the same time, the resonant spin- 
scattering allows electrons to pass through the SMM: 
the Kondo effect for transport 0, @ results in a zero 
bias conductance anomaly that has been studied experi- 
mentally in many systems with small spins (e.g. q uantum 
dots 0, M, IHIll 111 and single molecules Ojig). Such 
an effect is unexpected in SMMs because the S > 1/2 
underscreened Kondo effect is suppressed by the easy- 
axis anisotropy barrier which freezes the spin along the 
easy axis. However, we find that even a weak transverse 
anisotropy induces a pseudo-spin-1/2 Kondo effect. The 
corresponding Kondo temperature is experimentally ac- 
cessible due to a compensation by the large value of the 
physical spin S. We perform a scaling analysis |16| for 
the effective pseudo-spin-1/2 model and verify the results 
by a non-perturbative numerical renormalization group 
(NRG) calculation [H E3 for the full large-spin Hamil- 
tonian. 

Model. We consider SMMs which can be described by 
the following minimal model in the limit of strong tunnel- 



coupling to electron reservoirs H = Hm + : 



H M = -DSl -\Y. B ^ + (!) 



H K = JS • s + y^ e ka a\ a 
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Here S z is the projection of the molecule's spin on the 
easy axis, chosen as z-axis, S± — S x ± iS y , and s — 

Hkk' J2aa> a {a(l T )<™' a k'<T' is the local electron spin in 
the reservoir (t is the Pauli-matrix vector). As usual, 
the electronic states \ka) are a linear combination of 
states of both physical electrodes 0, H • The first term 
in Eq. Q describes the easy-axis magnetic anisotropy 
of the molecule, i.e. the states \S Z ), S z = —S, . . . , S are 
its eigenstates. The second term in Eq. Q describes 
transverse anisotropy perturbations which in general re- 
duce the symmetry to that of a discrete group of sym- 
metry operations caused by the geometrical structure of 
the molecule and its ligands via spin-orbit effects. The 
individual transverse terms written in our model are in- 
variant under 2n-fold rotation (n = 1, 2, 3) about the easy 
axis. To keep the notation systematic we deviate from 
the conventional notation for the anisotropies E = B2 
and C = B4. Note that the relative strength of the per- 
turbations is B 2n S 2{ - n - 1) /D. The first term in Eq. @ 
describes the exchange coupling of the molecular spin to 
the effective reservoir [the second term in Eq. and 
thereby transfers charge. The coupling is antiferromag- 
netic, J > 0, in the blockade regime as may be readily 
shown from the Schrieffer- Wolff transformation [l^ . 
We will show that in a half-integer spin SMM the Kondo 
effect lifts the blockade of both spin-tunneling (due to 
TR symmetry) and electron-tunneling (due to energy 
and charge quantization). Concerning the former effect, 
for half-integer S > 1/2 the eigenstates of the molec- 
ular Hamiltonian Hm are at least two-fold degenerate 
(Kramers doublets) and are linear combinations of states 
from only one of the disjoint sets {\^ z S±(2n)k)}k=o,i,2,..., 
see Fig. ^ Hence the transverse perturbations B^n can 
not connect the opposite magnetic basis states | ± S) as 
they do for integer spin S, i.e. QTM is blocked |20| . 
However, a Kondo spin-flip process, Eq. @, can change 
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FIG. 1: Coupling scheme of the ground states including QTM 
B 2 for 5 = 3/2,5/2,7/2. Black (dashed) and red (full) lines 
denote \ S Z ) states belonging to the different TR-invariant sub- 
spaces. 



S z by one and connect the disjoint sets of molecular 
states. Thus in cooperation with the QTM terms the 
molecular spin can be completely reversed. This is simi- 
lar to the co-tunneling of nuclear and electronic spins [2l| . 
Note that in a Kondo effect due to a positive easy- 
axis anisotropy (D < 0) was studied and no transverse 
anisotropics were considered. 

In the following we compare situations where either low- 
or high-symmetry QTM perturbations dominate. This 
can be achieved experimentally by a chemical modifica- 
tion of the ligands |23l | or in transport experiments by 
changing the binding of the molecule to the electrodes, 
which can be controlled mechanically in some setups J24[ . 

Poor-man scaling analysis. To describe the low energy 
properties of the above model we perform a poor-man 
scaling analysis 0|. This approach leads to similar re- 
sults as the full NRG calculations (shown in Fig. [210] and 
discussed at the end), but allows for a detailed discussion 
of the processes leading to Kondo physics. We truncate 
the spectrum to the twofold degenerate ground state |±) 
and obtain an effective spin-1/2 Kondo model 



Hctf = J 



\v){v'\{v\S\ 



V )s 



(3) 



vv' — ± 



with the pseudo-spin operators P± — P x ± iP y = |±)(t| 
and P z = (|+)(+| - |-)(-|)/2. The effective exchange 
constants depend on i?2n, D and <S* through 



j z = 2J(+\S Z \+) > 

j x ,y = J(+\S+±S-\-) 



(4) 
(5) 



which we have calculated numerically Importantly, these 
constants are completely anisotropic, except for special 
cases. The scaling equations are (2W — conduction elec- 
tron bandwidth, p — density of states) : 



dj a 
d\nW 



(6) 



where a, (3, 7 are cyclic permutations of x, y, z [23 . Spec- 
ification of any two scaling invariants j a 2 ~ji3 2 , ct 7^ (3, de- 
fines a 3-dimensional scaling curve. Inversion of any pair 
of j a ,jp leaves the scaling equations invariant, whereas 
inverting a single one reverses the flow. Interestingly, all 
scaling trajectories flow to the strong coupling limit ex- 
cept those in planes of uni-axial symmetry, \j a \ = < 



|j' 7 | with jajpj'y < 0. In the latter case one has a fer- 
romagnetic fixed line which is unstable with respect to 
infinitesimal perturbations perpendicular to it which are 
typically present in our model. If the effective exchange 
constants lie close to this line the Kondo temperature will 
thus be strongly suppressed. If the Kondo effect occurs 
and \j z \ > \j x \ > \j y \ with \j z \ ^ \j y \, we find for the 
Kondo temperature (defined here as the scale where the 
first coupling constant diverges) 



cs 



\n(T K /W cS ) = 



\jy\ 



3* -3x 

Jz Jy 



(7) 
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Here cs~ 1 (it|m) is the inverse of the elliptic integral 
cs(u|m), see |2g. In the uni-axial planes \j z \ = \j x \ 
or \jx\ = \j y \ Eq. (7J reduces to the well-known expres- 
sions for easy-axis anisotropy [53, 13 since cs _1 (it|0) = 
arctan(l/u) and cs~ 1 (u|l) = arctanh(Vl + u 2 ). In this 
model the upper bound of the Kondo scale Tk is the en- 
ergy separation between Kramers-degenerate ground and 
first excited state of the isolated molecule W e s(D, {B 2n }) 
which is of the order 0.1 meV ~ 1 K. 
According to Eq. J5J, the exchange couplings j X:V are gen- 
erated by spin-tunneling. This gives j z > \j x \, \j y \ for not 
too strong QTM. Thus, the only case where the Kondo 
effect can not be observed is \j x \ = \j y \ and j x j y < 0, 
which using Eq. (JSJ gives — ) = 0. This means 

that the Kondo effect is not observable when the spin 
raising operator of the original molecular spin can not 
flip the pseudo-spin from the down to the up value, an 
intuitively quite obvious condition. Most importantly, as 
we will illustrate in the following, this condition can be 
checked very easily provided the spin and the symmetry 
of the molecular magnet are given. If a P>i n quantum 
tunneling term is present, we get 



^0 
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(8) 



as a condition for the observability of the Kondo effect 
in molecular magnets with weak QTM. 
We first consider the limit of a dominant low-symmetry 
QTM term, B 2 > B 4 S 2 ,B 6 S 4 . We always find a spin- 
1/2 Kondo effect, see Fig.|2 because the three couplings 
are different except for B2 = D. At that point Eq. 
can be rewritten as Hm = 2DS y + const. The result- 
ing uni-axial symmetric couplings \j z \ — \j x \ > \j y \ al- 
low for a flow to the strong coupling fixed point, also 
in this case. The Kondo temperature, shown in Fig. [2 
has a non-monotonic dependence on B%/ D, which is en- 
hanced with increasing S. For B2 <C D (weak QTM), 
the criterion (JSJ applies and is always fulfilled for any 
half-integer spin. The states forming the two ground 
states |±) are connected by the spin raising operator 
by S — 1/2 QTM processes (each contributing a factor 
oc B2/D) and one co-tunneling process, see Fig.^ There- 
fore, the B2/D dependence, estimated from Eqs. 1)415(1 . is 
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FIG. 2: Kondo temperature, deduced from the NRG level 



flow with J = 0.1 W and D = 
QTM B 2 for 5 = 3/2, 5/2, 7/2. 



5 X 10 W, as function of 



\ Jz \ « 25 J » |j x | « 
gives for the Kondo temperature 



(B 2 /D) s -i. Using this 



(9) 



The exponent becomes 5- independent for S> 1. How- 
ever, the complicated spin-dependent prefactor left out 
in Eq. 10 decreases with 5 stronger than W e g increases: 
in this regime the Kondo temperature therefore decreases 
with increasing spin. Interestingly, this tendency changes 
for larger quantum tunneling. Near B 2 = D the perpen- 
dicular couplings dominate and grow with increasing 5: 
j XtX = Jy/S{S + 1) + 1/4 and j y = J for B 2 = D. From 
Eq. Q one obtains by expanding in j y /j z oc 1/5 an 
enhancement of Tr with S> 1: 



T, 



B 2 =D 
K 



/w t 



off ■ 



(10) 



In this expression there are two competing factors: in- 
creasing B 2 /D enhances the r.h.s. of Eq. J7J which is 
maximal at B 2 — D, but simultaneously reduces the 
splitting W e s between ground and excited state from 
(25 — 1)D to AD. Hence the maximal Kondo temper- 
ature occurs for a value B 2 /D < 1. Finally, for B 2 ^> D 
the Kondo temperature is suppressed again since one 
coupling becomes dominant: \j x \ » J ^> \j y \,\j z \- We 
note that in the case B 2 — D the molecule can not 
be considered as a molecular magnet and it is only dis- 
cussed here to explain the tendency of the increase of the 
Kondo temperature with increasing spin for strong quan- 
tum tunneling. We suggest to study SMMs with moder- 
ate quantum tunneling, as e.g. Feg 111 ' 1 with D — 0.27 K 
and B 2 = 0.046 K 29] . For the above mechanism to be 
relevant the spin 5 = 10 needs to be changed to an half- 
integer value by changing the charge via a gate electrode. 
The value of the Kondo exchange coupling depends on 
the details of the adjacent charge states, e.g. changes in 
anisotropics and total spin. For a quantitative calcula- 
tion further input from experiment and ab-initio calcu- 
lations (e.g. [30| ) is needed. 

Now we consider a dominant QTM perturbation of higher 
symmetry, B 4 > B 2 /S 2 ,B 6 S 2 . For B 2 = B 6 = 
we have 4 disjoint subsets of basis states which cannot 
be connected by the QTM term, e.g. for B 4 S 2 <C D 




FIG. 3: Scheme for the spin selection rule, (a) It is fulfilled 
for QTM B 4 and S = 5/2 leading to a Kondo effect, (b) For 
QTM B 4 and S = 7/2 the product of the effective couplings 
j x jy < and the Kondo effect is suppressed, (c) For QTM 
Be and spin S = 9/2. The Kondo coupling can not couple 
the ground state doublet. 



the ground states are linear combinations of { | =F S ± 
4fc)}fe=o,i,2,..., see Fig.^a,b). While condition © is ful- 
filled for S = 5/2 + 2m, (to = 0, 1, . . . ), Fig. ^a), it 
is violated for spin S = 7/2 + 2m. In the latter case, 
only the spin lowering operator can increase the pseudo- 
spin, Fig.|3Jb). However, with increasing B4S 2 /D a level 
crossing between the ground and excited state of different 
symmetry results in a sharp change in the Kondo tem- 
perature. Hence B4 induces a quantum-phase transition, 
see Fig. Hfor 5 = 7/2, c.f. 0. 

Finally, we mention the possibility that the Kondo cou- 
pling can not connect two TR-invariant subsets which 
leads to a complete vanishing of the effective coupling 
constants j XjV = 0. The lowest order QTM where this 
effect takes place is of sixth order: in general the Kondo 
coupling cannot overcome the "mismatch" |A5 2 | = 6, 
see Fig. [3Jc) . In this case only for intermediate or strong 
Bq a level crossing may give rise to ground states which 
support a Kondo effect. 

Numerical RG. We use Wilsons's NRG [T3II1 to check 
the results of our scaling analysis taking into account the 
full model (i.e. no truncation to a ground state dou- 
blet is made). As input parameters we take J = 0.1IT, 
number of states N = 1500, discretization A = 2 and 
D = 5 x 10 _3 IT. Since the original procedure has been 
formulated for a spin-1/2 Kondo model we have modi- 
fied it to incorporate an arbitrary spin of the impurity. 




FIG. 4: Kondo temperature, deduced from the NRG level 
flow, as function of QTM B 4 for S = 5/2, 7/2, 9/2. All other 
parameters as in Fig. 
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We analyzed the RG level flow as function of iteration 
number N- ltel in order to determine the low-temperature 
fixed point and the Kondo temperature which is defined 
as the energy scale where the crossover to strong cou- 
pling takes place. For half-integer spin S and D > we 
observe a flow to the strong coupling fixed point only 
for a QTM perturbation Bm 0, as expected from 
the above scaling analysis. The Kondo temperature for 
spin S = 3/2, 5/2, 7/2 and dominant B 2 QTM is plotted 
in Fig. |21 It shows in good qualitative agreement with 
the scaling results the discussed non-monotonic behav- 
ior as function of the QTM. For dominant B4 QTM the 
Kondo temperature is plotted for spin S = 5/2,7/2,9/2 
in Fig.0] The mapping onto a pseudo-spin- 1/2 system is 
valid as long as there is no crossing of levels and Tk does 
not exceed the gap to the first excited state. The former 
is the case for the experimentally most relevant regime 
B 2n S 2 ( n -V < D. 

Discussion. Consequently the observation of Kondo- 
tunneling through SMMs requires a judicious selection of 
three quantities: (i) the total spin must be half-integer; 
(ii) the dominant QTM perturbation should be moderate 
B2 n S 2n ~ 2 < D and (iii) the total spin should satisfy the 
selection rule © if a high-symmetry QTM term domi- 
nates. In experiments where the charge state, and hence 
the spin state, can be controlled by a gate voltage, the 
Kondo effect can only be seen in every (2n) subsequent 
Coulomb diamond for molecules with a dominant (2n) 



QTM perturbation. This assumes that the spin S in- 
creases/decreases by i every next charge state and has to 
be contrasted to the even/odd alternation of the Kondo 
effect usually observed in quantum dots. It would be 
of interest to mechanically alter the symmetry of the 
magnetic core of the SMM in situ, e.g. in a mechan- 
ically controlled break-junction setup |2J], and thereby 
suppress/enhance the Kondo effect. 

Summary. Using scaling and numerical renormaliza- 
tion group techniques we have found that spin- and 
electron-tunneling become correlated in half-integer spin 
magnetic molecules which are strongly coupled to elec- 
trodes. The spin-1/2 Kondo anomaly in the linear con- 
ductance signals the externally induced tunneling of the 
magnetization of the molecule. The Kondo tempera- 
ture shows a non-monotonic dependence on the relative 
strength of the transverse magnetic anisotropy of the 
molecule. Importantly, the large spin of SMMs is found 
to compensate for anisotropy effects which are expected 
to suppress Kondo physics. The symmetry of the trans- 
verse anisotropy imposes a selection rule: the Kondo ef- 
fect only occurs for selected values of the molecular spin. 
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